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Abstract 

We consider a hard jet production tagged by a muon pair in ultrarelativistic heavy 
ion collisions. The process cross section is calculated by the CompHEP Monte-Carlo 
generator taking into account full 7* jZ interference pattern at LHC energies. We have 
found that reasonable statistics, ~ 1000 events per 1 month of LHC run with lead beams, 
can be expected for realistic geometrical acceptance and kinematic cuts. The transverse 
momentum imbalance due to interactions of jet partons in the medium is evaluated for 
/i + /x~pair+jet correlation, as well as for the correlation between /i + \iT pair and a leading 
particle in a jet. Theoretical and experimental uncertainties of these observables are 
discussed. 
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1 Introduction 



One of the important tools to study properties of quark-gluon plasma (QGP) in ultrarelativistic 
heavy ion collisions is a QCD jet production. Medium-induced energy loss of energetic partons, 
the so-called jet quenching, has been proposed to be very different in cold nuclear matter and in 
QGP, resulting in many challenging observable phenomena [1]. Recent RHIC data on suppres- 
sion of inclusive high-p T charge and neutral hadron production from STAR [2], PHENIX [3], 
PHOBOS [4] and BRAHMS [5] are in agreement with the jet quenching hypothesis [6]. However 
direct event-by-event reconstruction of jets and their characteristics is not available in RHIC 
experiments at the moment, while the assumption that integrated yield of all high-pr particles 
originates only from jet fragmentation is not obvious. 

At LHC a new regime of heavy ion physics will be reached at ^/snn — 5.5 TeV where 
hard and semi-hard QCD multi-particle production can dominate over underlying soft events. 
The initial gluon densities in Pb— Pb reactions at LHC are expected to be significantly higher 
than at RHIC, implying stronger partonic energy loss which can be observable in various new 
channels [7, 8, 9]. In particular, the potentially important process is production of a single jet 
opposite to a gauge boson in 7+jet [10] and Z+jet [11] or a virtual photon in 7*(^ Z + Z~)+jet [12] 
final states, dominantly through processes such as 

q g ^ q 7 , qg^qZ, q g q 7* • 

In heavy ion collisions, the relative pr between the jet (or leading particle in a jet) and the 
boson becomes imbalanced due to interactions of the jet partons in the medium. 

In the 7+jet case the main problem arises from the jet pair production background when a 
leading n in the jet is misidentified as a photon. The "photon isolation" criteria usually used in 
pp collisions do not work with the same efficiency in high multiplicity heavy ion interactions [13]. 
Thus, the question of adequate using photon+jet correlation to study jet quenching requires 
further investigation. On the other hand, the production of jet tagged by dileptons is not 
affected significantly by backgrounds and can be used to observe p T -imbalance as a signal of 
medium-induced partonic energy loss. 

In this Letter we analyze dimuon+jet production (including both j*/Z — > modes) in 

heavy ion collisions at the LHC. In Sect. 2 the cross section of this process is calculated by 
CompHEP Monte-Carlo generator and the expected event rate is estimated for realistic geo- 
metrical acceptance and kinematic cuts. Sect. 3 describes shortly the model of partonic energy 
loss in QGP used to calculate p^-imbalance between n + fi~ pair and jet (or leading particle in a 
jet). Discussion on numerical results and their experimental and theoretical uncertainties can 
be found in Sect. 4, summary - in Sect. 5. 
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2 Dilepton+jet production at LHC 



We use the CompHEP Monte-Carlo generator [14] for cross section calculation of dilepton+jet 
process and subsequent generation of event for this process. CompHEP is a general, tree level 
generator, which allows one to study almost all processes 2 — > iV (up to N = 5) in the framework 
of the usual technique of Feynman diagrams squared for different models (SM, MSSM, many 
other...). It generates, squares and symbolically calculates a set of Feynman diagrams for a 
given process and creates a numeric Monte-Carlo generator for the process. This MC generator 
allows one to compute cross sections (with applied cuts), to build distributions and to generate 
events with partons in the final state; the initial partons are convoluted with parton distribution 
functions (PDF). 

For simplification of our calculation and event generation we apply a special hash-model in 
CompHEP [15]. In this model a unitary rotation of down quarks transfers Cabibbo-Kabayashi- 
Maskawa (CKM) matrix elements from interaction vertices to parton distribution functions. 
It allows one to unify two light quark generations to one only. As a result, this trick reduces 
significantly number of subprocesses which we need to take into account for the process. This 
model applies two approximations: 

• Uh (u, c) and dh (d, s) quarks are massless; 

• b and t quarks do not interact with light quarks and CKM matrix element V t b—1. 

In the problem of dilepton+jet investigation the influence of factors violating these conditions 
is very small so we can soundly use this approximation. 

There are nine subprocesses contributed to the process of the dilepton+jet production in 
the framework of SM with Zias/i-approximation. Feynman diagrams for the subprocesses are 
depicted on fig.l. We use the following physics parameter values: a = 1/127.9, sin6w = 0.2311, 
M z = 91.1876 GeV, T z = 2.4368 GeV, M M = 105.7 MeV, M b = 4.85 GeV, PDF are taken 
from cteq51 [16]. We apply a following set of "loose" cuts for event generation. In the further 
investigation these cuts have been strengthened: 

• Vt > 5 GeV/c and # > 20 GeV/c; 

• | rf> jet |< 3. 

We do not apply a cut on M(// + //~) because of a singularity in the region M(// + // - ) — > is 
regularized by massiveness of muon. Contributions of all subprocess to the total cross section 
is presented in the Table 1. 
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Table 1: Contributions to process pp — > ji fi + + jet (t/s^ = 5.5 TeV) 



Subprocess: 


Cross section (pb): 


u h u_ h -> + 5 


27.5 


44 -> + g 


18.1 


u h g -> + «/» 


105.0 


4# -> + 4 


36.2 




44.3 


d h g — > + 4 


21.9 




0.6 


6g — > + 6 


3.6 


6g — > + 6 


3.6 


gg -> n~n + + q 


260.9 



Note that the cross section for the process /i + yU~+jet was estimated in ref. [12] at RHIC and 
LHC energies taking into account the 7* contribution only. Such approximation seems valid for 
RHIC, but for LHC the contribution from Z and 7*/^ interference term are also significant. 

Now we estimate the expected event rate for realistic geometrical acceptance and kinematic 
cuts. To be specific, the geometry of Compact Muon Solenoid (CMS) detector is considered [17, 
18]: pseudo-rapidity coverage | i] \< 3 for jets and | rj \< 2.4 for muons. Extra cuts P^ +M > 50 
GeV/c and > 50 GeV were applied. Then the corresponding pp cross section for /i + /i~+jet 
production is w 16 ph, and Pb— Pb cross section is estimated as 16 pb x (207) 2 0.7 /ib. 
The corresponding event rate in a one month Pb— Pb run (assuming 15 days of data taking), 
T = 1.3 x 10 6 s, with designed luminosity L = 10 27 cm -2 s -1 , is iV cv = Ta PbPb L 1000 in this 
case. Note that potential using also e + e~+jet channel could increase observed event rates by a 
factor ~ 2 and requires further investigation. 

To conclude this section, let us discuss the potential background sources in the CMS ex- 
perimental situation. Semileptonic heavy quark decays and uncorrelated pion and kaon decays 
are expected to give main contributions to the dimuon spectra at LHC energies, ~ 10 5 events 
in a one month Pb— Pb run for CMS acceptance [8, 18]. However, the request of hard enough 
cut on muon pair transverse momentum, P^ +tJl > 50 GeV, together with the additional trigger 
to have a hard jet with > 50 GeV in the opposite hemisphere, makes such background 
"contamination" negligible. Moreover, the experimental control on background extraction may 
be done by monitoring uncorrelated and correlated sources independently. The uncorrelated 
part can be subtracted using like-sign dimuon mass spectra, while the correlated part can be 
rejected using tracker information on the dimuon vertex position [8, 18]. 
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3 Simulation of jet quenching at LHC 



In order to generate the initial distributions of jets and fi + fi~ pairs in nucleon-nucleon sub- 
collisions at y/s = 5.5 TeV, we have used CompHEP package for initial parton configuration 
setting and PYTHIA_6.2 [19] for subsequent jet fragmentation. After specifying initial partonic 
state, event-by-event Monte-Carlo simulation of rescattering and energy loss of partons in 
QGP is performed (for details of the model one can refer to [20, 21]). The approach relies 
on accumulative energy losses, when gluon radiation is associated with each scattering in the 
expanding medium together with including the interference effect by the modified radiation 
spectrum dE/dl as a function of decreasing temperature T. The basic kinetic integral equation 
for the energy loss AE as a function of initial energy E and path length L has the form 

v ' 

where I is the current transverse coordinate of a parton, dP/dl is the scattering probability 
density, dE/dl is the energy loss per unit length, A = \j{ap) is in-medium mean free path, 
p oc T 3 is the medium density at the temperature T, a is the integral cross section of parton 
interaction in the medium. Such a numerical simulation of the free path of a hard jet in QGP 
allows any kinematic characteristic distributions of jets in the final state to be obtained. 

The collisional energy loss due to elastic scattering with high-momentum transfer have orig- 
inally been estimated by Bjorken in [22], and recalculated later in [23] taking also into account 
the loss with low-momentum transfer dominated by the interactions with plasma collective 
modes. Since the latter process contributes to the total collisional loss without the large factor 
~ ln(E/fir)) (/Ad is the Debye screening mass) in comparison with high-momentum scatter- 
ing and it can be effectively "absorbed" by the redefinition of minimal momentum transfer 
t ~ n 2 D under the numerical estimates, we used the collisional part with high-momentum trans- 
fer only [21], 

3TE/2 
dE co1 1 r ' da 

and the dominant contribution to the differential cross section 

da_ ^ r 2na 2 s (t) 12tt 

dt ' as (33-2N f )\n(t/A 2 QCD ) {6) 

for scattering of a parton with energy E off the "thermal" partons with energy (or effective mass) 
m ~ 3T <C E. Here C = 9/4, 1,4/9 for gg, gq and qq scatterings respectively, a s is the QCD 
running coupling constant for Nf active quark flavors, and Aqcd is the QCD scale parameter 



which is of the order of the critical temperature, Aq C d — T c ~ 200 MeV. The integrated cross 
section a is regularized by the Debye screening mass squared jJ? D {T) ~ Arra s T 2 (l + Nj/6). 

There are several calculations of the inclusive energy distribution of medium-induced gluon 
radiation from Feyman multiple scattering diagrams. The relation between these approaches 
and their main parameters were discussed in details in the recent writeup of the working group 
"Jet Physics" for the CERN Yellow Report [7]. We restrict to ourself here by using BDMS 
formalism [24]. In the BDMS framework the strength of multiple scattering is characterized by 
the transport coefficient q = fi%/X g (X g is the gluon mean free path), which is related to the 
elastic scattering cross section a (3). In our simulations this strength in fact is regulated mainly 
by the initial temperature T . Then the energy spectrum of coherent medium-induced gluon 
radiation and the corresponding dominated part of the radiative energy loss has the form [24]: 



where T\—L/ (2A s ), y — uj/E is the fraction of the hard parton energy carried by the radiated 
gluon, and Cr = 4/3 is the quark color factor. A similar expression for the gluon jet can be 
obtained by substituting Cr = 3 and a proper change of the factor in the square bracket in (4), 
see ref. [24]. The integral (4) is carried out over all energies from uj min = E LPM = pL 2 D X g , the 
minimal radiated gluon energy in the coherent LPM regime, up to initial jet energy E. 

The medium was treated as a boost-invariant longitudinally expanding quark-gluon fluid, 
and partons as being produced on a hyper-surface of equal proper times r [25]. In order to 
simplify numerical calculations (and not to introduce new parameters) we omit the transverse 
expansion and viscosity of the fluid using the well-known scaling solution due to Bjorken [25] 
for a temperature and density of QGP at T > T c ~ 200 MeV: 



Let us remark that the influence of the transverse flow, as well as of the mixed phase at T = T c , 
on the intensity of jet rescattering (which is a strongly increasing function of T) seems to be 
inessential for high initial temperatures T 3> T c [20, 21]. On the contrary, the presence of 
viscosity slows down the cooling rate, which leads to a jet parton spending more time in the 
hottest regions of the medium. As a result the rescattering intensity goes up, i.e., in fact the 
effective temperature of the medium is increased as compared with the perfect QGP [20, 21]. 
We also do not take into account here the probability of jet rescattering in nuclear matter, 
because the intensity of this process and the corresponding contribution to the total energy 







e{t)t^ = e r 4/3 , T(r)r 1 /3 = T r 1/3 , p{r)r = p r . 



(6) 
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loss are negligible due to the much smaller energy density in "cold" nuclei. For certainty 
we used the initial conditions for the gluon-dominated plasma formation expected for central 
Pb-Pb collisions at LHC [26]: r ~ 0.1 fm/c, T ~ 1 GeV, p g « 1.95T 3 . For non-central 
collisions we suggest proportionality of the initial energy density e to the ratio of the nuclear 
overlap function and the effective transverse area of nuclear overlapping [21]. 

In each event the distribution of jet production vertex at the given impact parameter b of 
A A collision is generated according to the distribution [21] 

^ = tmm 

(tWCLV max 

J dip J rdrT A {r-i)T A (r- 2 ) 
o o 

where r 12 (b, r, ip) are the distances between the nucleus centers and the jet production vertex 
V(r cosip, rsin^); r max (b,ip) < Ra is the maximum possible transverse distance r from the 
nuclear collision axis to the V; Ra is the radius of the nucleus A; Ta(ti :2 ) is the nuclear 
thickness function (see ref. [21] for detailed nuclear geometry explanations). After that, in 
every ith scattering of the co-moving particle (with the same longitudinal rapidity) a fast 
parton loses energy in the collisions and radiatively, Aej = tj/(2m ) + uJi, where tj and ooi are 
simulated according to eqs. (2) and (4) respectively. Thus in each event the energy of an initial 
parton decreases by the value AE(r,ip) = J2i^ e i- 

In the frame of this model and using above QGP parameters we evaluate the mean energy 
loss of quark of E T = 50 GeV in minimum-bias Pb— Pb collisions, (AE^) ~ 25 GeV. In order 
to analyze the sensitivity of dimuon-jet correlations to the absolute value of partonic energy 
loss, we also performed the same calculations for the reduced initial temperature, T = 0.7 
GeV, which results in decreasing average energy loss by a factor ~ (1/0. 7) 3 ~ 3. 

The distribution over a difference between Pjt^ , transverse momentum of fi + fi~ pair, 
and Ej. , observed transverse energy of jet, depends crucially on a fraction of the partonic 
energy loss falling outside the jet cone. There are some discussions on the angular spectrum of 
in-medium radiated gluons in the literature [20, 24, 27, 28, 29]. In fact, since coherent Landau- 
Pomeranchuk-Migdal radiation induces a strong dependence of the radiative energy loss of a 
jet on the angular cone size, it will soften particle energy distributions inside the jet, increase 
the multiplicity of secondary particles, and to a lesser degree, affects the total jet energy. On 
the other hand, the energy loss in the collisions turns out to be practically independent of the 
jet cone size and causes the loss of the total jet energy, because the bulk of the "thermal" 
particles knocked out of the dense matter by the elastic scatterings flies away in on almost 
transverse direction relatively to the jet axis [20]. Thus, although the radiative energy loss of 
an energetic parton dominates over the loss in the collisions by up to an order of magnitude, 
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the relative contribution of the latter jet energy loss grows with increasing the jet cone size 
due to the essentially different angular structure of loss for two mechanisms [20]. Moreover, 
the total energy loss of a jet will be sensitive to the experimental capabilities to detect low-p T 
particles - products of soft gluon fragmentation: thresholds for a giving signal in calorimeters, 
influence of the strong magnetic field, etc. [18]. 

Since the full treatment of the angular spectrum of emitted gluons is rather sophisticated 
and model-dependent [20, 24, 27, 28, 29], we considered two simple parameterizations of the 
distribution of in-medium radiated gluons over the emission angle 6. The "small-angular" 
radiation spectrum was parameterized in the form 



where 9 ~ 5° is the typical angle of the coherent gluon radiation estimated in [20]. The 
"broad-angular" spectrum has the form 

dN 9 1 

We believe that such a simplified treatment here is enough to demonstrate the sensitivity of 
p T -imbalance in +jet production to the medium-induced partonic energy loss. 

4 Numerical results and discussion 

To be specific, the jet energy is defined here as the total transverse energy of the final particles 
collected around the direction of a leading particle inside the cone R = \/ Arj 2 + Aip 2 = 0.5, 
where r] and tp are the pseudorapidity and the azimuthal angle respectively. Fig. 2 shows the 
distribution over [P^ ^ — E^^j for the cases without and with medium-induced energy loss 
obtained in the framework of our model in minimum-bias Pb— Pb collisions, two parameteri- 
zations of distribution on gluon emission angles (8) and (9) being used. The same geometrical 
acceptance and kinematic cuts as in Sect. 2 were applied: | r? jct |< 3, | rf |< 2.4, pt^ > 5 
GeV/c, P^ +M ",^ et > 50 GeV. Although (p£ + ^ - E^) -distribution is already smeared in 
pp case, mainly due to the initial state gluon radiation, the mean value (^Pj> M — E^^, as 
well as the maximum of the distribution are close to 0. The partonic energy loss in heavy 
ion collisions results in the visible asymmetry of this distribution, its smearing and shifting 
mean and maximum values. The effect is more pronounced for the "broad-angular" radiation, 
because the contribution of the "out-of-cone" partonic energy loss is larger as compared with 
the "small-angular" radiation case. It is important to note that the contribution to jet-dimuon 
P T -imbalance from collisional part (always "out-of-cone" ) is rather dominant over the contribu- 
tions from "broad-angular" and "small-angular" radiation. The latter (dash-dotted histrogram) 




(8) 
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does not disappear totally just due to the fact that not only leading (parent) parton, but all 
partons of a jet pass through the dense medium and emit gluons under the angles 6 relatively to 
their proper directions, which in general may not coincide with the jet axis (determined by the 
direction of a leading particle) and sometimes be even at the jet periphery. Thus dimuon-jet 
correlation will be affected strongly by the collisional part of the loss. 

In the real experimental situation the jet observables will be sensitive to the accuracy of 
jet energy reconstruction in a high multiplicity environment. There are the two terms deter- 
mining this accuracy: systematic jet energy loss and jet energy resolution (due to calorimeter 
and jet finder peculiarities, influence of magnetic field, etc.). Short measure of jet energy will 
be well-controlled systematic error which has the similar values for heavy ion and for pp colli- 
sions, and it can be taken into account using the standard calibration procedure [30]. On the 
other hand, the finite jet energy resolution will result in additional smearing (Pj, ^ — E^\- 



distribution without shifting mean and maximum values. In order to illustrate the sensitivity 
of such observables to experimental jet energy resolution <te t , we set for definiteness the simple 
parameterization, oe t = 1.5y/Er GeV, which is close to one obtained with window-type jet 
finding algorithm for central Pb-Pb collisions (with background generated by fast MC code 
HYDRO [31] for multpilicity dN ± /dy(y = 0) = 8000) at CMS [7]. Fig.3 shows the same distri- 
bution as fig. 2, but final jet energy was smeared in each event by value oe t - Due to additional 
smearing the initial distribution, the effect of jet quenching on shifting mean and maximum 
values becomes less visible and rather marginal for moderate partonic energy loss. 

Since jet observables are affected by a number of theoretical (in particular, sensitivity to 
the angular spectrum of medium- induced radiation and to the collisional part of energy loss) 
and methodical (finite jet energy resolution) uncertainties, complementary leading particle 
measurements are potentially important, such as p^-imbalance between muon pair and a leading 



particle with maximum E T ) in a jet. The choice of a factor of 5 denotes only the fact that 
the most probable value of a fraction of jet energy carried by leading particles is ~ 0.2 in our 
case. This distribution is originally more smeared and asymmetric than the distribution over 
However additional smearing and shifting mean and maximum values due to 
partonic energy loss can be also clearly seen even for relatively small loss, (AE?) ~ 8 GeV 
at T = 0.7 (b = 0) GeV. Moreover, the observed p T -imbalance between pair and a 

leading particle in a jet is directly related to the absolute value of partonic energy loss and 
almost insensitive to the form of the angular spectrum of emitted gluons and collisional loss. 
The small difference for various angular distributions is just due to the moderate distinction 
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of event samples which were triggered by having a jet with > 50 GeV. Since fixing the 
minimum threshold for jet detection is sensitive to the total jet energy (and consequently 
to the angular dependence of energy loss), the small influence of the latter on dimuon- leader 
correlation appears. We also become convinced that taking into account the jet energy smearing 
(7e t almost does not have influence on dimuon-leader correlation (see fig. 5). 

5 Conclusions 

The channel with dimuon tagged jet production in ultrarelativistic heavy ion collisions was 
analyzed. The cross section of this process and corresponding event rates at LHC energies were 
evaluated with CompHEP Monte-Carlo generator taking into account full 7* jZ interference 
pattern. The reasonable statistics, ~ 1000 events per 1 month of LHC run with lead beams, 
can be expected for realistic geometrical acceptance and kinematic cuts. 

The correlations between pair and jet, as well as between pair and a leading 

particle in a jet, were first numerically studied for heavy ion collisions. The medium-induced 
partonic energy loss can result in significant smearing the distribution on difference between 
the transverse momentum of fi + fi" pair and the jet transverse energy, and shifting mean and 
maximum values of the distribution. This effect will be sensitive to the fraction of partonic 
energy loss (dependent on the form of the angular spectrum of in-medium radiated gluons) 
falling outside the jet cone. However the finite experimental jet energy resolution can result 
in additional smearing dimuon-jet correlation, which can make difficult observation of Pt- 
imbalance especially for moderate partonic energy loss. 

Since jet observables are affected by a number of theoretical and methodical uncertainties, 
complementary leading particle measurements will be useful and even preferable. P^-imbalance 
between fj, + fj,~ pair and a leading particle in a jet is quite visible even for moderate loss, directly 
related to the absolute value of partonic energy loss, almost insensitive to the angular spectrum 
of emitted gluons and to experimental jet energy resolution. 

Finally, the study of /i. + /i~+jet and /i + /z~-leading-particle correlations is important for ex- 
tracting the information about medium-induced partonic energy loss and properties of super- 
dense QCD matter to be created in heavy ion collisions at the LHC. 
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Figure 2: The distribution over difference between the transverse momentum of pair, 
Pj. ^ , and the jet transverse energy, Ej. , without (solid histogram) and with medium- induced 
partonic energy loss for the "small-angular" (8) (dotted histogram - radiaitve and collisional loss, 
dash-dotted historgam - radiaitve loss only) and the "broad-angular" (9) (dashed histogram) 
parameterizations of emitted gluon spectrum in minimum-bias Pb— Pb collisions. Applied 
kinematical cuts are described in the text. Initial QGP temperature T = 1 {b = 0) GeV (a) 
and T = 0.7 (b = 0) GeV (b). 
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Figure 4: The distribution over difference between the transverse momentum of pair, 
M , and five times transverse energy of a leading particle in a jet, 5 x _E^ adcr . The other 
conditions are the same as for fig.2. 
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